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The  Reynolds  number  similarity  hypothesis  of  Townsend  \The  Structure  of  Turbulent  Shear  Flow 
(Cambridge  University  Press,  Cambridge,  UK,  1976)]  states  that  the  turbulence  beyond  a  few 
roughness  heights  from  the  wall  is  independent  of  the  surface  condition.  The  underlying  assumption 
is  that  the  boundary  layer  thickness  S  is  large  compared  to  the  roughness  height  k.  This  hypothesis 
was  tested  experimentally  on  two  types  of  three-dimensional  rough  surfaces.  Boundary  layer 
measurements  were  made  on  flat  plates  covered  with  sand  grain  and  woven  mesh  roughness  in  a 
closed  return  water  tunnel  at  a  momentum  thickness  Reynolds  number  Re9  of  ~  14  000.  The 
boundary  layers  on  the  rough  walls  were  in  the  fully  rough  flow  regime  {kf  3s  100)  with  the  ratio  of 
the  boundary  layer  thickness  to  the  equivalent  sand  roughness  height  S/ks  greater  than  40.  The 
results  show  that  the  mean  velocity  profiles  for  rough  and  smooth  walls  collapse  well  in  velocity 
defect  form  in  the  overlap  and  outer  regions  of  the  boundary  layer.  The  Reynolds  stresses  for  the  two 
rough  surfaces  agree  well  throughout  most  of  the  boundary  layer  and  collapse  with  smooth  wall 
results  outside  of  3 ks.  Higher  moment  turbulence  statistics  and  quadrant  analysis  also  indicate  the 
differences  in  the  rough  wall  boundary  layers  are  confined  to  L<  5 ks.  The  present  results  provide 
support  for  Townsend’s  Reynolds  number  similarity  hypothesis  for  uniform  three-dimensional 
roughness  in  flows  where  8/kfs^  40.  [DOI:  10.1063/1.1843135] 


I.  INTRODUCTION 

The  rough  wall  boundary  layer  is  of  great  engineering 
interest,  including  the  applications  of  flow  over  ship  hulls 
and  airplanes,  as  well  as  fluid  transport.  Current  engineering 
models  treat  roughness  as  a  small  perturbation  to  the  smooth 
wall  boundary  layer.1  However,  if  the  effect  of  surface 
roughness  propagates  into  the  outer  layer,  the  application  of 
rough  wall  models  that  rely  on  smooth  wall  similarity  scaling 
will  yield  erroneous  results.  Understanding  the  extent  of  this 
perturbation  for  a  variety  or  roughness  types  would  improve 
modeling  and  predictive  capabilities. 

The  engineering  importance  of  predicting  the  effect  of 
wall  roughness  on  boundary  layer  and  pipe  flows  has  long 
been  identified.  Much  of  the  seminal  work  in  this  area  was 
carried  out  by  Nikuradse  and  Colebrook  and  White.  Later 
investigations  focused  on  the  effect  of  roughness  on  turbu¬ 
lence  structure. 4-6  An  extensive  review  of  the  present  knowl¬ 
edge  of  rough  wall  boundary  layers  is  given  by  Raupach  et 
al.  and  Jimenez. 

The  mean  velocity  profile  in  the  overlap  and  outer  re¬ 
gions  of  a  turbulent  boundary  layer  over  a  smooth  wall  is 
given  as 

U+=-ht(y+)  +  C+-wK),  (1) 

where  the  von  Karman  constant  k=  0.41,  the  smooth  wall  log 
law  intercept  C=5.0,  and  II  is  the  wake  parameter.  Both 
Clause/1  and  Hama10  pointed  out  that  the  primary  effect  of 
surface  roughness  on  the  mean  flow  is  to  cause  a  downward 
shift  in  this  profile.  For  “fc-type”  rough  walls,  the  downward 
shift  AU+  called  the  roughness  function,  correlates  with  k+, 


the  roughness  Reynolds  number.  The  mean  velocity  profile 
in  the  overlap  and  outer  regions  of  a  turbulent  boundary 
layer  over  a  rough  wall  can,  therefore,  be  stated  as 

U+  =  —  ln(y+)  +  C  -  AU+  +  —  w(  -  ) .  (2) 

K  K  \SJ 

By  evaluating  Eqs.  (1)  and  (2)  at  y=S,  Hama11  showed  that 
AU+  can  be  found  as  the  difference  between  the  smooth  wall 
log  law  intercept  C  and  the  rough  wall  intercept  C-AU+  at 
the  same  Re^.  The  roughness  function  is,  therefore,  found  as 


The  implicit  assumption  in  Eq.  (3)  is  that  the  mean  flow  for 
both  smooth  S  and  rough  R  walls  obeys  the  velocity  defect 
law  in  the  overlap  and  outer  regions  of  the  boundary  layer: 

u:-u+=f(^j.  (4) 

There  is  a  great  deal  of  experimental  support  for  a  universal 
velocity  defect  law,7'1 1-14  and  it  is  consistent  with  the  idea  of 
Reynolds  number  similarity  given  by  Townsend15  that  turbu¬ 
lence  outside  the  inner  layer  is  unaffected  by  surface  condi¬ 
tion.  Some  prominent  rough  wall  studies  claim  that  changes 
to  the  mean  flow  do  occur  in  the  outer  layer. 1617  This  is 
primarily  observed  as  an  increase  in  the  strength  of  the  wake 
II,  which  the  researchers  assert  is  due  to  the  higher  rate  of 
entrainment  for  rough  walls. 

There  is  also  a  large  degree  of  similarity  in  the  turbulent 
stresses  for  rough  and  smooth  walls.  The  Reynolds  number 
similarity  hypothesis  of  Townsend15  and  subsequent  exten- 
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sions  by  Perry  and  Chong  and  Raupach  et  al.  states  that 
the  turbulent  motions  are  independent  of  surface  condition 
outside  the  roughness  sublayer  at  sufficiently  high  Reynolds 
number.  This  implies  that  the  turbulent  stresses,  normalized 
by  the  wall  shear  stress,  are  universal  outside  of  the  rough¬ 
ness  sublayer.  The  dominant  view  at  present  is  that  the 
roughness  sublayer  extends  ~3k-5k  from  the  wall,  where  k 
is  the  roughness  height.  The  outer  region  has  been  generally 
thought  to  be  unaffected  by  the  roughness.  ’  This  is  sup¬ 
ported  by  a  number  of  studies,  including  the  work  of  Perry 
and  Li19  on  expanded  mesh  surfaces  and  Acharya  et  al.1"  for 
machined  surfaces.  Andreopoulos  and  Bradshaw  present 
both  mean  and  turbulence  profiles  for  flow  over  one  sandpa¬ 
per  surface.  Reynolds  stress  results  showed  good  collapse 
outside  the  roughness  sublayer.  However,  effects  in  the  triple 
products  were  noted  up  to  10k  from  the  wall.  Ligrani  and 
Moffat  ~  presented  mean  and  turbulence  results  for  closely 
packed  spheres  in  the  transitionally  rough  regime.  For  this 
type  of  roughness  element,  the  turbulence  quantities  col¬ 
lapsed  with  smooth  wall  results  outside  of  the  roughness 
sublayer.  Schultz  and  Flack13  examined  fully  rough  flow 
over  closely  packed  spheres  with  and  without  the  addition  of 
a  secondary  roughness  scale.  Quadrant  analysis  and  the  ve¬ 
locity  triple  products  indicated  that  the  changes  in  the  turbu¬ 
lence  structure  in  the  rough  wall  boundary  layers  were  con- 
fined  to  y<8ks.  Schultz  and  Flack  also  tested  two  sand 
grain  surfaces  of  varying  roughness  heights.  Their  results 
indicate  collapse  of  the  mean  velocity  in  defect  form  and 
collapse  of  the  normalized  Reynolds  stresses  profiles  for 
both  the  smooth  and  rough  surfaces  in  the  overlap  and  outer 
regions  of  the  boundary  layer. 

Other  researchers  have  observed  roughness  effects  well 
into  the  outer  layer.  The  studies  of  Krogstad  et  al . 16  for  wo- 
ven  mesh  roughness  and  Keirsbulck  et  al.  for  transverse 
bar  roughness  indicate  that  the  wake  strength  of  the  mean 
velocity  profile  is  increased.  Additionally,  Krogstad  and 
Antonia  noted  that  the  turbulence  structure  is  altered  in  the 
outer  layer  for  a  woven  mesh  roughness.  Antonia  and 
Krogstad14  also  observed  changes  in  the  Reynolds  shear 
stress  profiles  for  surfaces  covered  with  transverse  rods. 
These  studies  imply  the  effect  of  surface  roughness  may 
propagate  well  into  the  outer  region. 

An  explanation  for  the  disparate  findings  is  most  likely 
due  to  the  “strong”  roughness  used  in  the  investigations 
where  effects  were  observed  in  the  outer  layer.  A  strong 
roughness  here  is  defined  as  a  surface  whose  equivalent  sand 
roughness  height  ’  ks  is  a  significant  portion  of  the  inner 
layer  thickness.  This  view  is  supported  in  the  recent  review 
article  of  Jimenez8  which  asserts  that  surfaces  where  S/k 
=£40  may  exhibit  roughness  effects  well  into  the  outer  layer, 
since  most  of  the  log-law  region  may  be  destroyed  by  the 
presence  of  the  roughness.  The  size  of  the  woven  mesh, 
traverse  rods,  and  tranverse  bars  used  in  the  studies  of 
Krogstad  et  al.,16  Antonia  and  Krogstad,14  and  Keirsbulck  et 
al}1  correspond  to  S/ks~  15,  S/ks~  8,  and  S/ks~l,  respec¬ 
tively,  while  the  values  of  S/k  in  these  investigations  were 
—48,  —46,  and  —26,  respectively.  In  the  study  of  Schultz 
and  Flack13  on  uniform  packed  spheres,  in  which  differences 
in  the  turbulence  structure  were  observed  out  to  8 ks,  S/ks, 


and  S/k  were  —30.  The  use  of  ks  to  define  the  extent  of  the 
roughness  sublayer  instead  of  k  itself  was  proposed  by 
Schultz  and  Flack  because  it  provides  a  common  measure 
of  the  influence  of  the  roughness  on  the  mean  flow.  The 
difficulty  with  using  ks  for  predicting  roughness  effects  is 
that  it  cannot  be  determined  a  priori  from  a  physical  measure 
of  a  generic  roughness.  The  authors  believe  that  S/ks  is  a 
more  appropriate  parameter  than  S/k  in  identifying  the  rela¬ 
tive  strength  and  extent  of  the  roughness  influence  on  the 
turbulence  in  the  boundary  layer. 

In  the  present  investigation,  the  flow  over  two  three- 
dimensional  rough  surfaces,  sandgrain  and  woven  mesh,  will 
be  presented.  These  roughness  types  were  selected  because 
both  surfaces  have  been  previously  studied1”'16  with  contrast¬ 
ing  results  related  to  the  extent  of  roughness  effects  in  the 
boundary  layer.  The  height  of  the  roughness  was  chosen  to 
produce  fully  rough  turbulent  boundary  layers  while  main¬ 
taining  conditions  in  which  the  boundary  layer  thickness  is 
large  compared  to  the  equivalent  sand  roughness  height 
{S/ks^A 0).  The  mean  velocity,  Reynolds  stress  profiles, 
higher  order  moments,  as  well  as  quadrant  analysis  for  flow 
over  these  surfaces  will  be  compared  to  smooth  walls.  Re¬ 
sults  will  focus  on  the  extent  that  the  roughness  effects  pen¬ 
etrate  the  boundary  layer. 

II.  EXPERIMENTAL  FACILITIES  AND  METHODS 

The  present  experiments  were  conducted  in  the  closed 
circuit  water  tunnel  facility  at  the  United  States  Naval  Acad¬ 
emy  Hydromechanics  Laboratory.  The  test  section  is  40  cm 
X40  cm  in  cross  section  and  is  1.8  m  in  length,  with  a  tun¬ 
nel  velocity  range  of  0-6.0  m/s.  The  current  tests  were  run 
at  a  tunnel  speed  of  —3.8  m/s  (Ret=5.1  X  106).  Further  de¬ 
tails  of  the  water  tunnel  facility  are  given  in  Schultz  and 
Flack.1”'13  Three  surfaces  are  tested  in  this  study.  One  is  a 
smooth  cast  acrylic  surface.  The  other  two  are  rough  sur¬ 
faces;  one  covered  with  80  grit  wet/dry  sandpaper  and  the 
other  with  woven  mesh  (center  line  spacing=1.0  mm,  wire 
diameter=0.16  mm).  The  maximum  peak  to  trough  rough¬ 
ness  heights  k  are  0.69  mm  and  0.32  mm,  for  the  sandpaper 
and  mesh  surfaces,  respectively.  The  test  specimens  were 
inserted  into  a  flat  plate  test  fixture  mounted  horizontally  in 
the  tunnel.  The  test  fixture  is  the  same  as  that  used  by 
Schultz  and  Flack.1”'1'  The  forward  most  200  mm  of  the 
plate  is  covered  with  36  grit  sandpaper  to  trip  the  developing 
boundary  layer.  The  use  of  a  strip  of  roughness  was  shown 
by  Klebanoff  and  Diehl”  to  provide  effective  boundary  layer 
thickening  and  a  fairly  rapid  return  to  self-similarity.  The  test 
specimen  mounts  flush  into  the  test  fixture  and  its  forward 
edge  is  located  immediately  downstream  of  the  trip.  The  flap 
was  set  at  small  angles  to  make  minor  adjustments  to  the 
upstream  streamwise  pressure  gradient.  A  schematic  of  the 
test  water  tunnel  test  fixture  is  shown  in  Fig.  1. 

The  boundary  layer  profiles  presented  here  were  taken 
1.35  m  downstream  of  the  leading  edge  of  the  test  fixture. 
Profiles  taken  from  0.80  m  to  the  measurement  location  con¬ 
firmed  that  the  flow  had  reached  self-similarity.  The  trailing 
150  mm  of  the  flat  plate  fixture  is  a  movable  tail  flap.  This 
was  set  with  the  trailing  edge  up  at  —4°  in  the  present  ex- 
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FIG.  1.  Schematic  of  flat  plate  test  fixture. 


periments  to  prevent  separation  at  the  leading  edge  of  the 
plate.  The  physical  growth  of  the  boundary  layer  and  the 
inclined  tail  flap  created  a  slightly  favorable  pressure  gradi¬ 
ent  at  the  measurement  location.  The  acceleration  parameter 
K  was  — 2.0X10-8  and  did  not  vary  significantly  between 
the  test  surfaces. 

Velocity  measurements  were  made  using  a  TSI  FSA3500 
two-component,  fiber-optic  laser  Doppler  velocimeter 
(LDV).  The  LDV  used  a  four  beam  arrangement  and  was 
operated  in  backscatter  mode.  The  probe  volume  diameter 
was  —90  /xm,  and  its  length  was  —1.3  mm.  The  viscous 
length  scale  v/uT  varied  from  5  f/m  for  the  rough  walls  to 
7  /l m  for  the  smooth  wall.  The  diameter  of  the  probe  vol¬ 
ume,  therefore,  ranged  from  13  to  18  viscous  lengths  in  the 
present  study.  The  LDV  probe  was  mounted  on  a  Velmex 
three-axis  traverse  unit.  The  traverse  allowed  the  position  of 
the  probe  to  be  maintained  to  +10  /rm  in  all  directions.  In 
order  to  facilitate  two-component,  near  wall  measurements, 
the  probe  was  tilted  downwards  at  an  angle  of  4°  to  the 
horizontal  and  was  rotated  45°  about  its  axis.  Velocity  mea¬ 
surements  were  conducted  in  coincidence  mode  with  20  000 
random  samples  per  location.  Doppler  bursts  for  the  two 
channels  were  required  to  fall  within  a  50  /as  coincidence 
window  or  the  sample  was  rejected. 

In  this  study,  the  friction  velocity  uT  for  the  smooth  sur¬ 
face  was  found  using  the  Clauser  chart  method  with  log-law 
constants  /r=0.41  and  5  =  5.0.  For  the  rough  walls,  uT  was 
obtained  using  a  procedure  based  on  the  modified  Clauser 
chart  method  given  by  Perry  and  Li.19  Further  details  of  the 
procedure  are  given  in  Refs.  12  and  26.  For  all  the  test  sur¬ 
faces,  the  total  stress  method  was  also  used  to  verify  uT.  This 
method  assumes  a  constant  stress  region  equal  to  the  wall 
shear  stress  exists  in  the  overlap  and  inner  layer  of  the 
boundary  layer.  If  the  viscous  and  turbulent  stress  contribu¬ 
tions  are  added  together,  an  expression  for  uT  may  be  calcu¬ 
lated  as  the  following  evaluated  at  the  total  stress  plateau  in 
the  overlap  and  inner  layer: 

1/2 

(5) 


dU  — 
v —  —  u  v 


_  fy 


III.  UNCERTAINTY  ESTIMATES 

Precision  uncertainty  estimates  for  the  velocity  measure¬ 
ments  were  made  through  repeatability  tests  using  the  proce- 
dure  given  by  Moffat.  LDV  measurements  are  also  suscep¬ 
tible  to  a  variety  of  bias  errors  including  angle  bias,  velocity 
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bias,  and  velocity  gradient  bias,  as  detailed  by  Edwards. 
Fringe  bias  results  from  the  inability  to  sample  scattering 
particles  passing  through  the  measurement  volume  at  large 
angles  since  several  fringe  crossings  are  needed  to  validate  a 
measurement.  In  this  experiment,  the  fringe  bias  was  consid¬ 
ered  insignificant,  as  the  beams  were  shifted  well  above  a 
burst  frequency  representative  of  twice  the  freestream 
velocity.  Validation  bias  results  from  filtering  too  near  the 
signal  frequency  and  any  processor  biases.  In  general  these 
errors  are  difficult  to  estimate  and  vary  from  system  to  sys¬ 
tem.  No  corrections  were  made  to  account  for  validation 
bias.  Velocity  bias  results  from  the  greater  likelihood  of  high 
velocity  particles  moving  through  the  measurement  volume 
during  a  given  sampling  period.  The  present  measurements 
were  burst  transit  time  weighted  to  correct  for  velocity  bias, 
as  given  by  Buchhave  el  al.  Velocity  gradient  bias  is  due  to 
velocity  variation  across  the  measurement  volume.  The  cor¬ 
rection  scheme  of  Durst  et  al.30  was  used  to  correct  it' .  The 
corrections  to  the  mean  velocity  and  the  other  turbulence 
quantities  were  quite  small  and  therefore  neglected.  An  ad¬ 
ditional  bias  error  in  the  v'  measurements  of  —2%  was 
caused  by  introduction  of  the  w'  component  due  to  inclina¬ 
tion  of  the  LDV  probe.  Bias  estimates  were  combined  with 
the  precision  uncertainties  to  obtain  the  overall  uncertainties 
for  the  measured  quantities  at  95%  confidence.  These  were 
calculated  using  the  standard  methods  outlined  in  Coleman 
and  Steele.  The  resulting  overall  uncertainty  in  the  mean 
velocity  is  +1%.  For  the  turbulence  quantities,  u'2,  v  '2,  and 
u'v'  the  overall  uncertainties  are  +2%,  +4%,  and  +7%,  re¬ 
spectively.  The  uncertainty  in  itT  for  the  smooth  walls  using 
the  Clauser  chart  method  is  +3%,  and  the  uncertainty  in  uT 
for  the  rough  walls  using  the  modified  Clauser  chart  method 
was  +5%.  The  uncertainty  in  uT  using  the  total  stress  method 
is  +6%  for  both  the  smooth  and  rough  walls.  The  uncertain¬ 
ties  in  the  boundary  layer  thickness  <5,  displacement  thick¬ 
ness  5*,  and  momentum  thickness  6  are  +7%,  +4%,  and 
+5%,  respectively. 

IV.  RESULTS  AND  DISCUSSION 
A.  Mean  flow 

The  experimental  conditions  for  each  of  the  test  cases 
are  presented  in  Table  I.  Significant  increases  in  the  physical 
growth  of  the  boundary  layer  were  noted  on  both  rough  sur¬ 
faces  compared  to  the  smooth  wall.  The  sandpaper  showed 
increases  of  21%,  58%,  and  44%,  while  the  mesh  had  similar 
increases  of  18%,  57%,  and  45%  in  S,  S*,  and  0,  respec¬ 
tively.  As  stated  previously,  uT  was  determined  using  both  the 
total  stress  and  Clauser  chart  methods.  These  results  are  also 
given  in  Table  I.  The  agreement  between  the  methods  is 
within  5%  in  all  cases.  The  friction  velocity  used  for  the 
normalization  presented  in  this  paper  was  obtained  using  the 
Clauser  chart  method  due  to  the  slightly  lower  uncertainty. 
However,  it  should  be  noted  that  since  uT  obtained  using  the 
total  stress  method  was  slightly  less  for  both  the  smooth  and 
rough  walls,  normalization  based  on  this  method  would  not 
significantly  change  the  comparisons  presented  herein. 

Figure  2  shows  the  mean  velocity  profiles  for  all  three 
surfaces  plotted  in  inner  variables.  The  smooth  wall  log  law 
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TABLE  I.  Experimental  test  conditions. 


Surface 

Ue 

(m/s) 

Re„ 

(m/s) 

Total 

stress 

«T 

(m/s) 

Clauser 

A  U+ 

8 

(mm) 

<5* 

(mm) 

8 

(mm) 

5k/ S 

5  k„/S 

Smooth 

3.77 

10  220 

0.137 

0.141 

26.4 

3.20 

2.53 

Sandpaper 
k= 0.69  (mm) 
k+=  134 
(^=100 

3.77 

14  340 

0.180 

0.185 

7.7 

31.9 

5.05 

3.64 

0.11 

0.08 

Mesh 

k= 0.32  (mm) 
k+= 64 
*+=138 

3.81 

14  120 

0.187 

0.196 

8.5 

31.1 

5.03 

3.66 

0.05 

0.11 

is  shown  for  comparison.  Both  rough  surfaces  display  a  lin¬ 
ear  log  region  that  is  shifted  by  Af/+  below  the  smooth  pro¬ 
file  indicating  an  increased  momentum  deficit  on  these  sur¬ 
faces.  The  sandgrain  and  mesh  surface  produced  roughness 
functions  of  At/+=7.7  and  8.5,  respectively.  It  is  of  note  that 
the  roughness  function  is  higher  for  the  mesh  surface  at 
nominally  the  same  unit  Reynolds  number  even  though  the 
roughness  height  k  is  less  than  half  that  of  the  sandpaper. 
This  clearly  reinforces  that  the  roughness  height  alone  is  not 
a  good  indicator  of  roughness  function.  Proper  scaling  pa¬ 
rameters  for  roughness  must  also  account  for  differences  in 
parameters  such  as  surface  texture,  roughness  density, 
and  roughness  slope.34  It  was  demonstrated  by  Furuya  and 
Fujita  that  not  only  the  roughness  height  but  also  the  pitch 
to  diameter  ratio  is  an  important  parameter  in  determining 
Al/+  on  woven  mesh  roughness.  At  present,  however,  there  is 
no  adequate  means  of  predicting  the  effect  of  a  generic  rough 
surface  on  the  mean  flow  (i.e.,  Al/+)  from  measures  of  the 
roughness  profile.  It  seems  plausible  that  since  k  is  not  a 
reliable  indicator  (at  least  among  different  roughness  types) 
of  the  roughness  effect  on  the  mean  flow,  it  is  likely  not  a 
proper  length  scale  to  define  the  extent  of  the  roughness  ef¬ 
fect  on  the  turbulent  motions  (known  as  the  roughness  sub¬ 
layer)  either.  Alternatively,  the  equivalent  sand  roughness 


FIG.  2.  Mean  velocity  profiles  in  wall  coordinates  (overall  uncertainty  in 
U+:  smooth  wall,  +4%;  rough  walls,  ±7%). 


height  ks  provides  a  common  measure  of  the  influence  of  the 
roughness  on  the  mean  flow  and  is  related  to  \U+  via  the 
following,  ’  where  B,  the  log-law  intercept  for  uniform 
sand  grain,  is  equal  to  8.5: 

At/+  =  —  ln(fcj)  +  C-B.  (6) 

K 

The  estimated  extent  of  the  roughness  sublayer  [5k J  8)  is 
given  in  Table  I.  Values  of  5k/ S  are  also  presented  for  com¬ 
parison. 

The  mean  velocity  profiles  for  all  test  cases  are  pre¬ 
sented  in  defect  form  in  Fig.  3.  Also  shown  for  comparison 
are  the  results  of  the  smooth  wall  direct  numerical  simulation 
(DNS)  by  Spalart36  at  Ret,=  1410.  The  velocity  defect  pro¬ 
files  exhibit  excellent  collapse  in  the  overlap  and  outer  re¬ 
gions  of  the  boundary  layer  supporting  the  notion  of  a  uni¬ 
versal  defect  profile  for  rough  and  smooth  walls,  as  proposed 
by  Clauser9  and  Hama.10  Similar  results  were  also  observed 
by  Acharya  et  al.  for  mesh  and  machined  surface  roughness 
and  by  Schultz  and  Flack  for  a  variety  of  roughness 
types.  Krogstad  et  al.16  and  Keirsbulck  et  al.11  found  that  the 
defect  profiles  for  boundary  layers  over  woven  mesh  and 
transverse  bars,  respectively,  did  not  collapse  with  smooth 


FIG.  3.  Velocity  defect  profiles  [overall  uncertainty  in  ( Ue—U)luT :  smooth 
wall,  +5%;  rough  walls,  ±7%]. 
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FIG.  4.  Normalized  streamwise  Reynolds  normal  stress  profiles  (overall 
uncertainty  in  u'2/u2:  smooth  wall,  ±5%;  rough  walls,  +7%). 
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FIG.  6.  Normalized  Reynolds  shear  stress  profiles  (overall  uncertainty  in 
—u'v' lu2:  smooth  wall,  ±8%;  rough  walls,  ±10%). 


wall  profiles.  This  was  due  to  an  increase  in  II  for  the  rough 
walls,  which  presumably,  resulted  from  the  higher  rate  of 
entrainment  of  irrotational  flow  at  the  outer  edge  of  the 
boundary  layer. 16  The  present  authors  believe  the  differences 
in  results  can  be  attributed  to  the  strong  roughness  used  in 
the  studies  of  Krogstad  et  al. 16  and  Keirsbulck  et  al. 17  In 
both,  the  ratio  of  the  boundary  layer  thickness  to  the  equiva¬ 
lent  sand  roughness  height  S/ks  for  the  rough  wall  flows  was 
=£15. 


B.  Reynolds  stresses  and  quadrant  analysis 

The  normalized  Reynolds  stress  profiles  ( u'2/u2T ,  v ,2l ir„ 
and  -u'v' hq)  for  the  test  surfaces  are  presented  in  Figs. 
4-6.  Also  shown  are  the  results  of  the  smooth  wall  DNS  by 
Spalart36  at  Refl=  1410.  The  approximate  extent  of  the  rough¬ 
ness  sublayer  for  the  sand  grain  and  mesh  surfaces  (5 kj  8)  is 
given  for  reference.  The  streamwise  Reynolds  normal 
stresses  u'2lu2T  (Fig.  4)  for  the  smooth  and  rough  wall  show 
excellent  agreement  throughout  the  overlap  and  outer  regions 
of  the  boundary  layer.  Near  the  wall  (y/<5=£0.05),  u'2/u2T  is 
significantly  lower  on  the  rough  walls.  The  region  of  the 


FIG.  5.  Normalized  wall-normal  Reynolds  normal  stress  profiles  (overall 
uncertainty  in  v'2/u2T:  smooth  wall.  ±6%;  rough  walls,  ±8%). 


roughness  influence  corresponds  to  y^3ks  on  both  the  sand¬ 
paper  and  the  mesh  surfaces.  The  disappearance  of  the  near¬ 
wall  peak  in  u'2  was  shown  by  Ligrani  and  Moffat""  to  be 
indicative  of  a  boundary  layer  in  the  fully  rough  flow  regime. 

The  wall-normal  Reynolds  normal  stresses  v'2/u2  (Fig. 
5)  also  indicate  good  collapse  in  both  the  overlap  and  outer 
regions  of  the  boundary  layer.  This  result  is  in  agreement 

12  13  19  22 

with  a  number  of  studies.  ’  ’  In  the  present  study,  the 
agreement  in  v'2/u2  is  observed  for  y5?ks  for  both  rough 
surfaces.  The  studies  of  Krogstad  et  al . 16  and  Keirsbulck  et 
al . 17  both  showed  that  significant  increases  in  the  wall- 
normal  Reynolds  normal  stress  penetrate  well  into  the  outer 
layer  over  rough  walls.  It  should  be  noted  that  even  though 
the  present  rough  surfaces  are  in  the  fully  rough  flow  regime, 
they  are  much  “weaker”  in  that  S/ks= 63  for  the  sandpaper 
and  S/ks= 45  for  the  mesh,  as  compared  to  S/ks =£  15  in  those 
studies. 

The  normalized  Reynolds  shear  stress  {-u'v' lit2)  pro¬ 
files  are  presented  in  Fig.  6.  Good  collapse  of  the  profiles  is 
again  observed  across  almost  the  entire  boundary  layer  (y 
^  ks) .  In  order  to  further  investigate  possible  differences  in 
the  flow  dynamics  between  smooth  and  rough  wall  boundary 
layers,  quadrant  analysis  was  carried  out  using  the  hyper- 
bolic  hole  size  H  method  of  Lu  and  Willmarth.  This  tech¬ 
nique  allows  the  contributions  of  ejection  Q2  and  sweep  QA 
motions  to  the  Reynolds  shear  stress  to  be  calculated.  The 
contribution  to  u'v'  from  a  given  quadrant  Q  can  be  ex¬ 
pressed  as 

_  1  fT 

(u'v')q=  hm  —  I  u'v'(t)IQ(t)dt,  (7) 

T-+  oo  l  Jo 

where  IQ{t )  is  a  trigger  function  defined  as 

|l  when  \u'v'\q  3=  H\lu'2yv'2 

Iq=  I  (8) 

[  0  otherwise . 

Figure  7  shows  the  normalized  contribution  from  ejection 
and  sweep  events  to  the  Reynolds  shear  stress  for  H=  0.  The 
profiles  of  both  the  Q2  and  QA  contributions  for  the  smooth 
and  rough  walls  show  excellent  agreement  across  almost  the 
entire  boundary  layer.  All  the  surfaces  display  a  nearly  linear 
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FIG.  7.  Normalized  Reynolds  shear  stress  contributions  from  Q2  and  Q4 
with  H= 0  (overall  uncertainty  in  —(u'v')m/u2  and  —{u'v')qJu\:±  10%). 
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FIG.  9.  Ratio  of  the  Reynolds  shear  stress  contributions  from  Q2  and  04 
with  H—2  (overall  uncertainty  in  (u' v!)q2I {u' v')q/]  :  ±  11%). 


decrease  in  the  Q2  and  Q4  Reynolds  stresses  in  the  outer 
layer.  This  indicates  that  roughness-induced  changes  to  co¬ 
herent  turbulent  structures  are  confined  to  a  near-wall  rough¬ 
ness  sublayer.  It  is  of  note  that  ejection  events  dominate 
sweeps  throughout  the  boundary  layer  on  the  smooth  wall 
and  for  y  3=  3ks  on  the  rough  walls.  The  results  of  Krogstad  et 
al.Ul  indicate  an  increase  in  both  the  Q2  and  Q4  contribu¬ 
tions  for  rough  walls  that  extends  well  into  the  outer  layer. 
Schultz  and  Flack  "  also  observed  an  increase  in  these  quan¬ 
tities  for  uniform  sphere  roughness,  but  only  for  y^5ks.  It 
should  be  noted  that  both  of  these  studies  were  carried  out  on 
stronger  roughness  (S/ks^  30)  than  the  present  investigation. 

Figure  8  shows  contribution  to  the  Reynolds  shear  stress 
from  strong  (//=  2)  ejection  and  sweep  events.  This  H  corre¬ 
sponds  to  instantaneous  Reynolds  shear  stress  events  larger 
than  5  u'v' .  The  profiles  of  both  the  Q2  and  Q  4  contributions 
for  the  smooth  and  rough  walls  again  show  excellent  agree¬ 
ment  outside  the  near-wall  region.  One  difference  that  can  be 
discerned  is  an  increase  in  strong  Q2  events  for  y/  0.05 
on  the  smooth  wall  that  is  absent  on  the  rough  walls.  Also, 
strong  Q4  events  decrease  on  the  smooth  wall  over  the  same 
range.  This  is  illustrated  more  clearly  in  Fig.  9,  which  shows 
the  ratio  of  the  contributions  from  Q2  and  Q 4  for  H=  2.  Near 
the  wall,  the  ratio  is  greater  than  unity  for  the  rough  walls 
and  less  than  unity  for  the  smooth  wall.  This  indicates  that 


0.15 


0.10 


'3 


0.05 


0.0  0.2  0.4  0.6  0.8  1.0  1.2  1.4 
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while  strong  ejection  events  play  a  larger  role  near  a  smooth 
wall,  sweeps  have  a  larger  contribution  near  rough  walls. 
This  was  also  observed  by  Krogstad  et  al. 16  and  is  likely  due 
to  the  difference  in  the  boundary  condition  at  y  =  0  for 
smooth  and  rough  walls.  On  a  smooth  wall,  v '  =0  at  y=0.  On 
a  rough  wall,  y=0  is  located  somewhere  between  the  rough¬ 
ness  peaks  and  troughs,  not  on  the  wall,  so  the  vertical  ve¬ 
locity  fluctuations  may  be  nonzero.  The  ratios  of  the  Q2  to 
Q4  contributions  presented  in  Fig.  9  show  good  agreement 
for  all  of  the  surfaces  for  y>3ks. 

C.  Velocity  triple  products  and  higher-order  moments 

While  the  present  Reynolds  stress  results  indicate  there 
is  a  great  deal  of  similarity  in  the  turbulence  on  smooth  and 
rough  walls  for  y>3ks,  higher-order  turbulence  moments 

9  j 

will  now  be  addressed.  Andreopoulos  and  Bradshaw” 
showed  that  the  velocity  triple  products  provide  a  sensitive 
indicator  of  changes  in  turbulence  structure  due  to  wall  con¬ 
dition,  however,  relatively  few  rough  wall  studies  have  pre¬ 
sented  these  statistics  with  the  exception  of  Antonia  and 
Krogstad,14  Kiersbulck  et  al.,11  Andreopoulos  and 
Bradshaw,21  and  Bandyopadhyay  and  Watson.39  The  distri- 


FIG.  8.  Normalized  Reynolds  shear  stress  contributions  from  Q2  and  Q4  FIG.  10.  Normalized  velocity  triple  product  u'O ul  (overall  uncertainty  in 

with  H=  2  (overall  uncertainty  in  —(u'v')q2Iu2t  and  —(u'v'Iq^u2:  ±12%).  u'2lu\:  ±28%). 
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y/S  y/5 


FIG.  11.  Normalized  velocity  triple  product  v 1  V/T  (overall  uncertainty  in 
V5lu\-.±  38%). 


FIG.  13.  Normalized  velocity  triple  product  u'v'2/u3  (overall  uncertainty  in 
u'v'2!u\ :  ±18%). 


butions  of  the  normalized  velocity  triple  products  u'3/u3T  and 
u'3lu\  are  presented  in  Figs.  10  and  11,  respectively.  The 
profiles  of  u,3/u3T  show  good  agreement  between  the  smooth 
and  rough  walls  for  y  >  5 ks.  Closer  to  the  wall,  some  differ¬ 
ences  can  be  noted.  The  sign  of  u,3/u3T  is  negative  for  the 
smooth  wall  for  y<0.05c>,  whereas,  it  is  positive  for  the 
rough  walls.  This  may  be  due  to  reduced  sweep  events  ( u ' 
>  0)  for  the  smooth  wall  than  for  the  rough  wall,  which  was 
also  identified  in  the  quadrant  analysis.  The  profiles  of 
v'3 / u\  also  show  agreement  within  their  experimental  uncer¬ 
tainty  outside  the  near-wall  region.  This  contrasts  with  the 
results  of  Antonia  and  Krogstad14  who  found  large  changes 
in  this  triple  product  for  flows  over  transverse  rod  roughness. 
This  was  observed  as  a  negative  value  of  vl3ln3T  over  a  sig¬ 
nificant  portion  of  the  boundary  layer  indicating  transport  of 
turbulent  kinetic  energy  towards  the  wall  instead  of  away 
from  it,  as  seen  for  a  smooth  wall. 

The  distributions  of  the  velocity  triple  products  ul2v' /u3T 
and  u'v'2/u3t  are  presented  in  Figs.  12  and  13,  respectively. 
The  quantity  u'2v'/u3T  represents  the  wall-normal  turbulent 
transport  of  u'2.  The  profiles  for  the  smooth  and  rough  walls 
agree  for  y>5ks.  Much  closer  to  the  wall,  there  is  a  large 


reduction  in  u'2v' lu3T  for  the  rough  walls.  A  similar  trend 
was  observed  for  mesh  roughness  by  Krogstad  and 
Antonia.40  The  gradient  du'2v' / dy  represents  the  turbulent 
diffusion  of  u'2  in  the  Reynolds  stress  transport  equation  for 
u’2.  On  the  smooth  wall,  there  is  a  loss  in  u'2  through  tur¬ 
bulent  diffusion  for  0.01  <y/ S< 0.025,  while  there  is  a  gain 
for  the  rough  walls.  The  normalized  turbulent  flux  of  Rey¬ 
nolds  shear  stress  u'v'2/u3T  also  indicates  good  agreement  for 
the  smooth  and  rough  walls  over  most  of  the  boundary  layer 
(Fig.  13).  The  differences  in  u'u'2/u3T  are  confined  to  the 
near-wall  region,  where  the  wall-normal  turbulent  transport 
of  Reynolds  shear  stress  (- u'v ')  is  toward  the  wall  on  the 
rough  walls  and  away  from  it  on  the  smooth  wall.  This  dif¬ 
ference  in  the  near-wall  transport  was  also  observed  by  An- 
dreopoulos  and  Bradshaw”  and  is  likely  due  to  the  stronger 
sweep  events  that  occur  on  the  rough  wall,  as  discussed  pre¬ 
viously. 

The  skewness  factor  distributions  for  u'  and  v' ,  Su  and 
Sv,  are  shown  in  Figs.  14  and  15,  respectively.  Also  shown 
for  comparison  are  the  smooth  wall  data  of  Antonia  and 
Krogstad14  at  Re0=  12  570.  There  is  excellent  collapse  of  ,S'„ 
for  all  three  surfaces  for  y  >  5 ks  and  good  agreement  with  the 
smooth  wall  data  of  Antonia  and  Krogstad14  outside  the 
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FIG.  12.  Normalized  velocity  triple  product  u'2v'/u3  (overall  uncertainty  in 
u'2v' lu3:  ±22%). 
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FIG.  14.  Skewness  factor  S„  (overall  uncertainty  in  Slt :  ±25%). 
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FIG.  15.  Skewness  factor  Sv  (overall  uncertainty  in  S„:  ±35%). 

roughness  sublayer.  Closer  to  the  wall,  the  rough  walls  show 
positive  skewness,  while  the  smooth  wall  has  negative  skew¬ 
ness.  The  positive  skewness  near  the  rough  walls  may  be  due 
to  the  less  strict,  wall-normal  boundary  condition,  allowing 
for  more  high  momentum  fluid  to  be  swept  into  the  near  wall 
region.  Sv  shows  agreement  within  the  experimental  uncer¬ 
tainty  across  the  entire  boundary  layer  (Fig.  15)  and  agrees 
well  with  the  smooth  wall  data  of  Antonia  and  Krogstad14  for 
y>5ks.  Antonia  and  Krogstad14  noted  significant  differences 
in  Sv  for  transverse  rod  roughness  well  into  the  outer  region 
of  the  boundary  layer.  The  flatness  factor  distributions  for  u' 
and  v',  Fu  and  Fv,  are  shown  in  Figs.  16  and  17,  respectively. 
Also  included  for  comparison  are  the  results  of  Bandyo- 
padhyay  and  Watson  for  flow  over  transverse  bars.  Both 
flatness  factors  show  good  collapse  for  smooth  and  rough 
surfaces  throughout  the  entire  boundary  layer  and  agree  rea¬ 
sonably  well  with  the  results  of  Bandyopadhyay  and  Watson 
over  most  of  the  boundary  layer. 

V.  CONCLUSIONS 

Comparisons  of  structure  of  turbulent  boundary  layers 
developing  over  two  roughness  types  and  a  smooth  wall  have 
been  made.  The  rough  surfaces  included  sandpaper  and  wo- 


FIG.  16.  Flatness  factor  Fu  (overall  uncertainty  in  F„:  ±3%). 


FIG.  17.  Flatness  factor  Fv  (overall  uncertainty  in  F„:  ±6%). 


ven  mesh;  surfaces  that  have  previously  been  used  to  exam¬ 
ine  the  concept  of  turbulence  similarity  in  the  outer  layer.12'16 
The  present  results  indicate  that  the  mean  velocity  profiles 
for  rough  and  smooth  walls  collapse  well  in  velocity  defect 
form  in  the  overlap  and  outer  regions  of  the  boundary  layer. 
The  Reynolds  stresses  and  quadrant  analysis  results  for  the 
two  rough  surfaces  agree  well  throughout  most  of  the  bound¬ 
ary  layer  and  collapse  with  smooth  wall  results  for  y>3ks. 
The  velocity  triple  products  and  higher  moment  turbulence 
statistics  also  indicate  that  the  differences  in  the  rough  wall 
boundary  layers  are  confined  to  y<  5 ks.  These  results  pro¬ 
vide  compelling  support  for  Townsend’s  Reynolds  number 
similarity  hypothesis  for  uniform  three-dimensional  rough¬ 
ness  in  flows  where  the  roughness  is  a  relatively  small  per¬ 
turbation  to  the  smooth  wall  case.  A  survey  of  the  literature 
along  with  the  present  results  indicates  that  S/ks  (not  S/k  as 
Jimenez8  recently  suggested)  is  the  proper  parameter  to  indi¬ 
cate  if  the  roughness  effect  on  the  turbulence  will  be  strong 
or  weak.  For  roughness  where  8/  ks  3=  40,  significant  similar¬ 
ity  in  the  turbulence  structure  can  be  expected  outside  the 
roughness  sublayer.  In  cases  where  SI ks  <  40,  turbulence 
modifications  may  be  anticipated  to  extend  well  into  the 
outer  layer.  Flows  with  strong  roughness  have  great  practical 
importance.  These  range  from  engineering  applications  of 
flows  through  fouled  heat  exchanger  tubes  and  over  fouled 
gas  turbine  blades  to  atmospheric  boundary  layers  in  urban 
or  forested  areas.  Therefore,  while  the  concept  of  turbulence 
similarity  appears  applicable  to  many  rough  wall  flows,  fur¬ 
ther  understanding  and  predictive  capability  are  needed  in 
cases  where  roughness  effects  are  large  and  turbulence  simi¬ 
larity  fails. 
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